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We study measurement-based quantum computation (MQC) using as quantum resource the planar 
code state on a two-dimensional square lattice (planar analogue of the toric code). It is shown that 
MQC with the planar code state can be efficiently simulated on a classical computer if at each step 
of MQC the sets of measured and unmeasured qubits correspond to connected subsets of the lattice. 



I. INTRODUCTION 

Quantum mechanical systems allow, for a class of com- 
putational problems, exponentially more efficient pro- 
cessing of information than classical systems. The ques- 
tion of where this speedup comes from has been under 
intensive debate over the recent years. 'Largeness of 
Hilbert space' pj, 'entanglement' (see e.g. @) and 'su- 
perposition and interference' for example, have been 
suggested. They all seem to be a part of the puzzle but 
it is difficult to pin down a single characteristic property. 

A prerequisite for a quantum speed-up in computation 
with a quantum system is the hardness of its classical 
simulation. One may learn about the cause for the quan- 
tum speed-up by investigating the circumstances under 
which it vanishes, i.e. when efficient classical simula- 
tion becomes available. A number of such scenarios have 
been described in the literature. For example, the uni- 
tary evolution of a spin chain can be followed efficiently as 
long as the entanglement with respect to all possible bi- 
partitions remains small Q- Further, linear optics with 
Gaussian states, systems of non-interacting fermions and 
qubit registers acted upon by gates from the so-called 
Clifford group are quantum systems which can be effi- 
ciently simulated classically; see 0, H, El and Q, re- 
spectively. 

Here we address classical simulation of quantum sys- 
tems in the context of measurement-based quantum com- 
putation (MQC), in particular the one-way quantum 
computer (QCc) &]■ In this scheme, one-qubit measure- 
ments are performed on a multi-qubit entangled resource 
state, the so-called cluster state. After the universal clus- 
ter state is created, no further interaction among the 
qubits takes place. Quantum information is written onto 
the cluster, processed and read-out from the cluster by 
the one-qubit measurements alone. 

To obtain a better understanding about MQC one may 
apply certain alternations to the original scheme and 
search for properties which remain invariant under these 
changes. Which other schemes for the processing of the 
measurement outcomes exist? Which other quantum 
states are universal resources and which properties char- 
acterize them? With regard to the classical processing of 



measurement outcomes in MQC, modified schemes have 
been described in 0, 0, |n| . Alternative universal re- 
sources for MQC have been presented in fill IT^L IT^| . For 
example, universal resource states exist in which each 
qubit is arbitrarily close to a pure state |ll| • Concerning 
universal vs. efficiently simulatable quantum resource 
states, a systematic study has begun in [l3L IbH IT5I Il6| | . 
Diverging amounts of entanglement, measured in terms 
of the so-called entanglement width, are required for uni- 
versality [T3L On the other side of the spectrum, 
MQC can be simulated efficiently classically by a so- 
called tree tensor network (TTN) if the resource state 
is a graph state and the graph is a tree or close to a tree 
|l4 llol [lfj . The prototypical example is the graph state 
on a line graph, i.e., the ID cluster state |l7|. An exam- 
ple for large deviation from tree-ness is the universal 2D 
cluster state whose TTN simulation is thus hard. 

In this paper, we describe a complementary simulation 
method for MQC, centered around planarity of graphs. 
The counterpart of the ID cluster state is the planar 
code state 113, and that of the TTN is the partition 
function of the Ising model. Within our framework, an 
example for large deviation from planarity again is the 
2D cluster state. 

Planar code states and cluster states are closely re- 
lated. For example, if one applies a certain pattern of 
Pauli-measurements to the two-dimensional cluster state, 
one can prepare the planar code state j2£j. One dimen- 
sion higher up, the fault-tolerance properties found in 
three-dimensional cluster states |2fJ are related to the 
'Random plaquette Z 2 gauge model in three dimensions' 
which also describes fault-tolerant data storage with a 
planar code [2l| . 

A further interesting property of the planar code state 
is that it obeys the entanglement area law j2^. That is, 
the entanglement entropy of a block of spins is propor- 
tional to its perimeter. Thus bi-partite entanglement in 
the planar code state is large. This state also exhibits 
topological quantum order and it is therefore not possi- 
ble to prepare the planar code state from a product state 
by a small-depth unitary quantum circuit |23| . 

However, our result is that MQC with the planar code 
state as the quantum resource is not universal and can be 
simulated efficiently classically. This unexpected prop- 
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FIG. 1: Stabilizer generators for the planar code state \K). 
Qubits are located on the edges of the lattice. Z-type stabi- 
lizer generators live on vertices, and X-type stabilizer gener- 
ators live on plaquettes. 



erty of the planar code state can be attributed to the 
exact solvability of the Ising model on a planar graph. 
Thus, although large entanglement in the resource state 
is necessary for MQC, it is not sufficient. 



II. THE PLANAR CODE STATE AND MQC 

We consider two-dimensional square lattice of dimen- 
sions L x L. It consists of N2 = L 2 plaquettes, Ni = 
2L(L + 1) edges, and N = (L + l) 2 vertices. Qubits 
live on the edges of the lattice, so the Hilbert space is 
H = (C 2 )®^ 1 . Let X e and Z e be the Pauli operators a x 
and a z acting on the qubit of edge e tensored with the 
identity operators on all other edges. For any vertex s 
and any plaquette p define stabilizer generators 



A, 



e£<5s 



B„ 



(i) 



Here 5s denotes the set of edges incident to vertex s, and 
dp denotes the set of edges making up the boundary of 
plaquette p, see Figure ^ Any pair of generators com- 
mute with each other since the sets Ss and dp always 
share even number of edges. All N2 plaquette-type gen- 
erators B p are independent. As for the vertex- type gen- 
erators A s , only No — 1 of them are independent, since 
FL=i ^ s = I- Thus the total number of independent 
generators is Nq + N2 — 1 = 2L(L +1). It coincides with 
the number of qubits Ni , so there exists a unique planar 
code state \K) G TL satisfying the stabilizer equations 



A S \K) = \K), B P \K) = \K), 



(2) 



for all s = I,..., No and p = 1, . . . , N2. It is worth- 
while to write down the state \K) explicitly. We shall 
use the computational basis of TL. Accordingly, a basis 
vector I a;) represents a configuration of N\ binary vari- 
ables x = {x e }e=i,...,AT 1 , x e € {0,1} living at the edges 
of the lattice. We shall call such a configuration a 1- 



even for every vertex s. In other words, a 1-cycle must 
have even number of edges incident to any vertex. Let 
Z C {0, l}^ 1 be a set of all 1-cycles. Then the state \K) 
can be explicitly written in the computational basis as 



\K) = 



1 



(3) 



One can easily check that \Z\ — 2 N2 , that is each plaque- 
tte contributes one independent 1-cycle to Z. 

MQC with the planar code state is defined as adaptive 

sequence of one-qubit measurements M]. 

applied to the state \K). After all measurements are 
done, every qubit on the lattice is measured in some basis. 
The choice of the measurement Mj + i may be determined 
by the outcomes of all earlier measurements Mi , . . . , Mj . 
We allow only complete von Neumann measurements. In 
other words, Mj is specified by a unity decomposition 

i=\^){^\ + \^))^)\ 

for some orthonormal basis \il>j), IV}) G C 2 . Let ej denote 
the qubit to be measured at step j and nij £ {0, 1} be the 
outcome of the measurement Mj . After the measurement 
the qubit ej is projected onto a state I^J* 3 )- Thus MQC 
is specified by 

• Description of the first measurement 
Mr = ( ei) IV?), W)) 

• Efficient algorithm that takes as input outcomes 
mi, . . . , nij of the first j measurements and returns 
a description of the next measurement 



M j+ i = (c i+ i,|^ +1 ) 



chain. Let us call a 1-chain x a I- cycle iff 



X e IS 



The output of MQC is a string of outcomes mi , . . . , m^ 
having some probability distribution p(mi, . . . , mjvj). 
We say that MQC can be efficiently simulated classically 
if there exists a classical randomized algorithm running 
in a time poly(L) that allows one to sample mi, . . . , m^ 
from the probability distribution p(rai, . . . , m jy 1 ) • 

In general, there are no restrictions on the order in 
which the qubits are measured in MQC. We shall con- 
sider more restrictive settings (because our technique 
does not apply to the general case). Consider a set 
Ej = (ei, e%, . . . , ej) including all qubits measured up 
to the step j of MQC. Let Ej be the complimentary set 
including all unmeasured qubits. Let us say that a set 
of edges E is connected iff for any pair of edges e, e' € E 
there exists a path eo = e, e\, e%, . . . , e„ = e' that con- 
nects e with e' and such that e; € E for all I. Our main 
result is 

Theorem. Suppose at each step j = 1, . . . , Ni of MQC 

the sets of measured and unmeasured qubits Ej , Ej are 
connected. Then MQC can be efficiently simulated clas- 
sically. 
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The proof relies on two observations. The hrst obser- 
vation made by Kitaev [24[ is that the overlap problem 
— computing the inner product between \K) and an ar- 
bitrary product state can be reduced to computing par- 
tition function of the Ising model with complex weights. 
The connection between overlaps of states and partition 
functions of models in statistical mechanics is explored 
in greater generality in 25] . 

More strictly, we use Barahona's techniques |26j to 
express the overlap as a generating function of perfect 
matchings on a properly defined planar graph. It was 
shown by Kasteleyn j27| that such generating functions 
can be computed efficiently, see for more recent ex- 
position. Thus if every qubit on the lattice is measured 
in some specified basis, probability of any particular out- 
come can be computed efficiently 

However, this observation solves only part of the prob- 
lem. First of all, since the number of possible outcomes, 
2 Nl , is exponentially large, computing probability of any 
particular outcome does not allow us to sample the out- 
comes according to these probabilities. More impor- 
tantly, since the measurements in MQC are adaptive, we 
must be able to compute probabilities of partial mea- 
surements, when only some subset of qubits is measured. 
Our second observation is that the overlap problem for a 
subset of qubits can be reduced to computing partition 
function of the Ising model for two independent replicas 
of the lattice with properly identified boundaries. The 
constraint that the subsets Ej , Ej are connected is the 
simplest way to ensure that the doubled lattice corre- 
sponds to a planar graph. The partition function of Ising 
model on any planar graph can be computed efficiently 
using Barahona's algorithm 26] . Given probabilities of 
any partial measurement outcomes, we can compute con- 
ditional probabilities and thus we can simulate MQC. 

Whether or not the connectivity constraint is a severe 
restriction on MQC is an open question. Note that this 
constraint is trivially satisfied in the standard scheme of 
MQC, when the qubits are measured in the order "from 
the left to the right". We also conjecture that the the- 
orem above can be generalized to the following scenar- 
ios: (1) the number of connected components in Ej,Ej 
is greater than one, but can be bounded by a constant. 
(2) The square lattice on a plane can be replaced by a 
lattice on any two-dimensional orientable surface with 
genus g bounded by a constant. The intuition support- 
ing these conjectures comes from the work psf providing 
an efficient algorithm to compute the generating func- 
tion of perfect matching if the genus of the graph can be 
bounded by a constant. 



III. CLASSICAL SIMULATION 

This section is organized as follows. Section UlI Al dc- 
scribes the reduction from the overlap problem to the 



Ising model and then to the perfect matchings on a pla- 
nar graph. It mainly follows the reference |26| (except 
for considering complex weights) and aims to make our 
exposition self-sufficient. Section Till Bl shows how to de- 
scribe a mixed state obtained from \K) by the partial 
trace over a subset of qubits as a mixture of planar code 
states. Finally, Section IIII CI explains how to compute 
probability of any partial measurement outcome by tak- 
ing two copies of the lattice. 

A. The overlap between planar code states and 
product states 

In this section we considerplanar code states on sub- 
graphs of the square lattice |3J] . We show that the over- 
lap (the inner product) between the planar code state 
and a product of any one-qubit states can be efficiently 
computed using techniques developed by Kasteleyn and 
Barahona pll^. 

Suppose the planar code state \K) is defined on a 
square lattice £, see Figure ^ Denote V{C) and E{C) 
sets of vertices and edges of C. For any subset of edges 
E C E{C) define a planar graph Ge = (V,E), where 
V C V(C) is a set of all vertices s £ V(C) having at least 
one incident edge from E, 

V = {se V(C) : dsn E ^ 0}. 

We shall denote C\(E) a set of 1-chains on the subset E. 

Let qubits live on edges e € E. Accordingly, the 
Hilbert space is now (C 2 )®I B I, and basis vectors corre- 
spond to 1-chains x £ C\{E). We shall consider a set of 
1-cycles 

Z(E) = {x £ Ci(E) : J2 x e [s even for a11 s e V }- 

eeSs HE 

(4) 

For brevity, we shall keep using notation Z = Z(E) 
throughout this section. Define a planar code state \Ge) 
associated with the graph Ge as 

i g *> = ^Ei*>- ( 5 ) 

VKI xez 

Note that \Z\ can be computed efficiently since Z is spec- 
ified by mod-2 linear equations. 

Consider an arbitrary product state 

l*> = ®l&>, \4>e) =a e \O)+0 e \l). 
eeE 

The quantity we are interested in is the overlap between 
|$) and \G E ), 

r=< G ,i*)=i*i-*n«.£ n (£)• 
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This expression becomes singular if a e = for some e. 
The singularity can be avoided by assigning small non- 
zero value to a e and using continuity of T as a function 
of a's and /3's (see also |35|). For every edge e £ E define 
a complex weight w e = e /a e . Then T is proportional to 
a "partition function" 

z =n n w - ( 6 ) 

x£lZ eGE:x e — l 

Let us mention that Z can be identified with a parti- 
tion function of the Ising model with complex weights. 
Indeed, introduce virtual Ising spins cry = ±1 living on 
faces / of Ge and use the ansatz 1 — 2x e = <Jfa g , where / 
and g are the two faces of Ge whose boundary includes 
the edge e. Then any choice of variables 07 automat- 
ically satisfies the 1-cycle constraint x G Z. Besides, 
w e" — \A^e exp (— Pf i9 (JfUg), where e^- 9 = yfw~e~. Thus 
Z is proportional to the partition function of the Ising 
models with Ising spins living on faces of Ge- In fact, 
the reverse reduction from the Ising model to the model 
described by Eq. is the first step in Barahona's so- 
lution of the Ising model |2^|. Since some arguments 
in Barahona's original work assume that the weights w e 
are real and positive, we shall briefly explain how to ap- 
ply the same approach to complex weights. (Note that 
in general partition functions with complex weights are 
more difficult to compute because of the sign problem.) 

Since Ge is a subgraph of the square lattice, vertices 
of Ge may have degree (number of incident edges) 1, 2, 3, 
or 4. Let us firstly get rid of vertices with degree 1, 2, 4. 
Suppose s £V has degree 1 and e is the edge incident to 
s. Clearly, x e = for any 1-cycle x £ Z. Thus Z does not 
depend on w e and we can safely remove the edge e from 
the graph. Now we can assume that all vertices of Ge 
have degree 2,3, or 4. Suppose s £ V has degree 2 and 
ex, e2 are the two edges incident to s. Clearly, x ei — x e2 
for any 1-cycle x £ Z. Therefore we can remove the 
vertex s and merge the edges ei, e2 into a single edge e 
with a weight w e = w ei w e2 . Now we can assume that all 
vertices of Ge have deg ree 3 or 4. Suppose s £ V has 
degree 4. Following 26] let us replace s by two vertices 
si, S2 of degree 3 by adding a new edge e = (s\, S2) with a 
weight w e — 1, see Figure|21 If ei , e2, e^, are the edges 
incident to s, see Figure [21 an d x £ Z is a 1-cycle, then 
x ei +x e2 +x e3 +x ei — (we use mod-2 arithmetics). After 
adding a new edge e we get x e = x ei + x ei = x e2 + x e3 . 
Thus there is one-to-one correspondence between 1-cycles 
on the original and the modified graph. Since w e = 1, the 
partition functions Eq. © are the same for both graphs. 
Now we can assume that all vertices of Ge have degree 
3. Let G = (V, E) be the 3-valent planar graph that 
we obtained by performing all the reductions above (by 
abuse of notations, we use the same symbols V ,E for 
new set of vertices and edges). We want to compute the 
partition function Eq. © for the graph G. Following |2^| 
let us introduce a new planar graph G = (V, E) such 




FIG. 2: Local modification of the graph proposed by Bara- 
hona [2^| to eliminate degree-4 vertices. The two graphs co- 
incide outside the dashed circle. The added edge e carries 
trivial weight w e = 1. The partition functions of the two 
graphs coincide. 




FIG. 3: Local modification of the graph proposed by Bara- 
hona |2ll to reduce summation over 1-cycles to summation 
over perfect matchings. The two graphs coincide outside the 
dashed circle. All six added edges e^,...,eg carry trivial 
weight, w e4 = . . . — w?g = 1. Any 1-cycle x £ Z on the 
graph G can be uniquely extended to a 1-chain y on G such 
that x and y coincide outside the dashed circle (including the 
edges ei,e2,e3) and such that y is a perfect matching inside 
the dashed circle. Conversely, y can not be a perfect matching 
if 2/ei + 2/e 2 + y e3 is odd. 

that each vertex s £ V with incident edges ei,e2, £ 
E is replaced by 6 vertices s\, . . . , s§ £ V and 6 edges 
64, . . . , eg as shown on Figure^ We assign trivial weights 
to all new edges, w ei = . . . = w eg = 1. A 1-chain y £ 
C\ (E) is a perfect matching on the graph G iff any vertex 
s £ V has exactly one incident edge e with y e = 1. Let 
M be a set of all perfect matchings on G. As was pointed 
out by Barahona [2(| , there is one-to-one correspondence 
Z(E) = M. between 1-cycles on G and perfect matching 
on G, in particular 

z=y: n ^=t, n ( ? ) 

x£Z e£E:x e = l yEM e eE:y^ = l 

By construction, G has at least one perfect matching 
(the one corresponding to the zero 1-cycle in G), and 
thus G must have even number of vertices, | | = 2m. 

Let A be the adjacency matrix of the graph G, i.e., 
A Stt — 1 if (s,t) £ E and A Stt = otherwise. By def- 
inition, A is symmetric, A T = A. One can rewrite the 
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summation over perfect matchings in Eq. J7J as 



2 m m! 

a-es 2m 3=1 

where we label vertices of G by integers 1,2,..., 2m, 
where a is a permutation of the 2m vertices, 



1 2 3 4 ■ • • 2m - 1 2m 

<7l er 2 <J 3 cr 4 ••• CT 2m _i (72m 



and S2m is the symmetri c gr oup. 

Kasteleyn has shown [23 that for any planar graph 
it is possible to find antisymmetric matrix B such that 
B s .t = iA s .t for any pair of vertices s, t, B T = —B, and 
such that B obeys the Pfaffian orientation condition: 

m m 
n^-i,^- =Sgn(cr) Y\_ ^0-23-1, 0-23- 
3=1 3=1 

for any permutation a € S2m- Here sgn(cr) is the parity 
of er. Define a complex antisymmetric 2m x 2m matrix 
B w such that Bf t = B s . t w (Sit) if (s, t) £ E and B™ t = 
otherwise. Then the partition function Eq. © can be 
expressed as the Pfaffian of B w , 

1 m 



j=l 



Now it can be computed efficiently using the well-known 
identity Pf(X) 2 = det (X) valid for any complex anti- 
symmetric matrix X. 



B. Tracing out qubits from the planar code state 

Suppose at some step j of MQC a subset of qubits 
E = {ex, . . . C i?(£) has been already measured. 
We want to show that a mixed state 

p E = Tr^ E \K)(K\ 

describing qubits from E only (before the measurement) 
can be represented as a probabilistic mixture of pure 
states simply related to the planar code state \Ge), see 
Eq. ©. 

Indeed, let E — E(C)\E be the set of traced out 
qubits. Denote V C V{C) a set of vertices having at 
least one incident edge from E, and dE = VD V be a set 
of vertices having at least one incident edge from both 
sets E, E (it can be regarded as a boundary of the set E). 
A 0-chain on some set of vertices M is a configuration of 
binary variables assigned to vertices of M , u = {u s }seMj 
u s 6 {0, 1}. Let us denote Cq(M) a set of all 0-chains 
on M. Given any 1-chain x € C\(E), define a bound- 
ary dx £ Co(V) as a 0-chain such that (dx) s is equal to 



the parity (0 or 1) of edges e incident to s with x e = 1. 
Define a subspace of relative 1- cycles 

Z(E, dE) = {x£ Ci(E) : (dx) s = for all s <£ dE}. 

Thus a relative 1-cycle satisfies the 1-cycle condition ev- 
erywhere except (may be) of vertices s £ dE. Define a 
subspace of syndromes as 

S(E) = {u£ Co(dE) : u = dz for some z £ Z(E, dE)}. 

Simple linear algebra arguments show that any relative 1- 
cycle z £ Z(E, dE) with a syndrome u = dz £ S(E) can 
be represented as z = z(u) + x, where z(u) £ Z(E,dE) 
is some fixed relative 1-cycle satisfying dz(u) — u and 
x £ Z{E) is 1-cycle. Similarly, one can define a subspace 
of syndroms 

S{E) = {u£ C (dE) : u = dz for some z £ Z(E, dE)}. 

Note that in general S(E) 7^ S(E). Any relative 1-cycle 
z £ Z(E, dE) with a syndrome u = dz £ S(E) can be 
represented as z — z(u) + y, where z(u) £ Z(E,dE) 
is some fixed relative 1-cycle satisfying dz(u) = u and 
y £ Z(E) is a 1-cycle. 

The above definitions and arguments imply that a 
Schmidt decomposition of \K) with respect to the parti- 
tion E(C) = E U E can be chosen as follows: 



\K) = ^= J2\Ge(u))®\G e (u)), 



\S 



(9) 



uES 



where 



S = S(E)nS(E), 



V\ Z ( E )\ xeZ (E) 



and 



vl^C-^Jl yez(E) 

Therefore we arrive at 

P E = W\ E \Ge(u))(G e (u)\. 



(10) 



C. Calculating probabilities of partial 
measurements 

In this section we show how to compute the overlap 
($|/0b|$), where |$) is an arbitrary product state. At 
this point we shall exploit the constraint that E and E 
are connected sets. 
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The first simplification arising from this constraint is 
that the subspace of syndromes S includes all even 0- 
chains on the boundary dE: 



S = {ueC (dE) 



Z^ ' 

s£dE 



= 0(mod2)}. (il) 



Indeed, consider arbitrary even 0-chain u G Co(dE), and 
let si,S2, ■ ■ ■ ,S2k € dE be the set of vertices such that 
u Sj = 1. Consider firstly the pair of vertices s\,S2- 
Since si,S2 € dE, there exist edges ei,e2 € -E inci- 
dent to Si,S2 respectively (may be e\ = €2)- Since E 
is a connected set, there exists a path 712 C E con- 
necting ei and e2- By definition, this path is a relative 
f-cycle, 712 € -Z(E, <9E) with the boundary c?7i2 sup- 
ported on si and S2. Similarly, using the connectivity 
of E, one can show that there exists a relative 1-cycle 
712 € Z(E,dE) with the boundary $712 supported on 
sx and S2- Repeating this arguments for the remaining 
pairs of vertices (S3, S4), . . . , (s2fc-ij S2/c) we can construct 
relative 1-cycles z(u) = 712 + ■ • • + 72fc-i,2fc € -Z(E, dE) 
and z(u) = 712 + . . . + 72fc— i,2fe € Z(E,dE) such that 
dz(u) = u and dz(u) = u. Thus ueS. Conversely, for 
any relative 1-cycle z 6 Z(E, <9E) the boundary dz must 
be even since z consists of elementary edges and every 
edge has two endpoints. The same is true for relative 
1-cycles on E. Therefore we have proved Eq. (|11|) . It 
follows that \S\ =2l aE l- 1 . 

The next simplification that we draw from the connec- 
tivity of E, E is 

Lemma. Suppose E, E are connected sets. Then the 
graph Ge can be drawn without intersections on a disk 
such that all vertices of dE lie on the boundary of the 
disk. 

Proof: Since Ge is a planar graph, any point of the 
plane either belongs to Ge, or is 'internal' with respect to 
Ge, or is 'external' with respect to Ge- Denote Ext (Ge) 
a set of all external points. The fact that E, E are con- 
nected implies that (1) The boundary dE is contained 
in the boundary of Ext (Ge)', (2) The region Ext (Ge) is 
topologically equivalent to a plane with a hole. There- 
fore one can smoothly deform a plane such that after the 
deformation is contained in a disk and all points of 
dE lie on the boundary of the disk. □ 

This lemma provides us with a way to compute the 
overlap ($|/3e|$). Indeed, consider two copies of the 
graph Ge drawn on disjoint disks D and D* . Let us 
denote the two copies as and G* E . Let us glue the 
disks D, D* into a sphere such that vertices of the bound- 
ary dE on the disk D are glued with the corresponding 
vertices of the boundary dE on the disk D* . Denote the 
resulting graph Ge^G* e . Since Ge^G* e can be drawn on 
a sphere, it is a planar graph and thus we can efficiently 
compute its partition function, as defined in Eq. ©. Let 
w e be a weight assigned to edge e e E on disk D. Let 



us assign the complex conjugated weight w* to the corre- 
sponding edge e on the disk D* . For any syndrome iie5 
let z(u) E Z(E,dE) be some fixed relative 1-cycle such 
that dz(u) = u. Then the partition function for Ge^G* e 
can be written as 



Z(G E UG E ) = J2 l Z («)! 5 



ues 



where 



Z(u)= J2 



n 



xGZ(E) eeE:(z(u)+x) c = l 

One can easily see that Z(u) is equal to the overlap be- 
tween a state \G E (u)) in the decomposition Eq. I|1(J|) 
and a product state |$) = ® e eE + w e that is 
Z(u) = y/\Z(E)\ (G E (u)\$). Therefore we conclude that 



(*|ps|*> = 



1 



2l» B l- 1 \Z(E)\ 



Z(G E U G* e ) 



and thus the overlap can be computed efficiently using 
techniques from the previous section. 

We can use the algorithm above for computing the 
overlap between pe and a product state to compute 
probabilities of any particular partial measurement out- 
comes p(m\, . . . , rrij). It allows us to compute condi- 
tional probabilities p(mj || mi, . . . , rrij—i) for any j = 
1, . . . , N\. If we already know the outcomes mi, . . . , rrij, 
we can sample rrij from the probability distribution 
p(m,j || mi, . . . , rrij-x). Thus all quantum steps in MQC 
can be efficiently simulated on a classical computer (with 
random numbers). 



IV. EXTENSION TO NON-PLANAR GRAPHS 

In this section we outline an MQC simulation scheme 
built around the planarity of grap hs. It is complemen- 
tary to the simulation scheme [T^ - 0| via tensor net- 
works, centered around the tree-ness of graphs. We iden- 
tify a parameter n which measures deviation from pla- 
narity and show that simulation of MQC is efficient in 
the number of particles in the resource state but expo- 
nentially inefficient in 77. 

In Section IIV Al it is shown that simulation based on 
planarity can not be subsumed under the tree tensor net- 
work method. In Section IIVBI we show that there is 
another interesting example among the states for which 
MQC can be related to the Ising model partition func- 
tion: a universal 2D cluster state. The corresponding 
interaction graph (to be distinguished from the graph 
upon which the definition of graph states is based) is 
highly non-planar, such that simulation — as expected — 
is not efficient. In Section HV CI we characterize the com- 
plexity of classical MQC simulation for non-planar Ising 
interaction graphs. 



7 



We would like to remark that the connection between 
the ID cluster state and non-interacting fcrmions, a topic 
interwoven with the solvability of the Ising model on pla- 
nar graphs, has previously been discussed in |29j . 

A. The planar code state and tensor contraction 

Here we show that a planar code state on a lattice of 
size L x L has an entanglement width of at least L such 
that its simulation scheme based on a TTN is exponen- 
tially inefficient in L. See [25| for an alternative proof of 
this result. 

The entanglement measure entanglement width has 
been defined in 0] and related to universality of MQC 
and to the complexity of its TTN-simulation in 0, ^| . 
We briefly restate the definition here. Consider a quan- 
tum system on a set Q of qubits. Be T a tree graph with 
each vertex having f or 3 incident edges (subcubic tree). 
The vertices with one edge are called leaves, and each 
of them corresponds to a qubit of the considered quan- 
tum system. If an edge e is deleted from T, the resulting 
graph T\e has two connected components which induce 
a bi-partition (A^,,B^) of Q. For a pure state we 
denote by ^A" ,B e (|^)) the entanglement entropy with 
respect to that bi-partition. The entanglement width is 
now defined as 

^wd(IV')) := minmaxE A e B e (lip)), 

T e 11 

where the minimization is over all subcubic trees with n 
leaves, and the maximum over all edges e in a given tree. 

Lemma. The entanglement width of a planar code state 
\Gi) on a square lattice of size L x L is 

E wd (\G L )) > L. (12) 

Proof: First we show that for any given subcubic 
tree T with \Q\ leaves, there exists an edge e such that 
1/3 \Q\ < \A%,\,\B%,\ < 2/3 \Q\. Given T, pick an arbi- 
trary 3-valent vertex as root. This induces an ancestor- 
descendant relation among pairs of vertices. For each 
vertex a £ V(T), denote by n(a) the number of leaves 
among a and all its descendants. There are vertices 
a € V(T) for which n(a) > 2/3n. Among those, there 
must be a vertex ao such that for both its direct descen- 
dants ai,a2 holds n(aq), n(a?) < 2/3 n. Choose a\ s.th. 
n(ai) > n(ci2). The desired edge then is e = (ao,ai). 
Be Aj. the set of qubits associated with the descen- 
dant leaves of a\. Then, \Aj\ — n{a\) < 2/3 n. Also, 
1/3 n < n(o )/2 = (rc(ai) + n(a 2 ))/2 < n(ai) = \A^\. 

Now, the entanglement entropy Ea b ,B e is linear in the 
length of the boundary between Aj. and [22^. It is 
minimized if A^ and are simply connected and use 
the external boundary of the L x L lattice. This occurs 
e.g. for Afp filling up the bottom part of the lattice up 




FIG. 4: Universal graph state and its surface code represen- 
tation, a) The two-dimensional cluster state \Cl), for L — 3. 
Qubits are located on the vertices », 0. b) Local unitary 
equivalent surface code state \Gl)- Here, qubits are located 
on the edges, and state stabilizers are associated with the sites 
and the plaquettes. To better illustrate the identification of 
qubits in \Gl) with their counterparts in \Cl), the edges of 
Gl are labeled by the same symbols as the vertices in Cl- 
The dashed lines indicate pairs of corresponding qubits. 

to a certain height. The entanglement entropy in that 
case is L or L + 1. Thus, max e Ea^,,b^, QGl)) > L, for 
all subcubic trees T. □ 
The planar code state considered here is, like any other 
stabilizer state, local Clifford equivalent to a graph state 
[3fl l3ll H3- For graph states, the entanglement width 
equals the so-called rank width 33] of the underlying 
graph. The rank width \ is the critical parameter for the 
complexity of MQC simulation via tree tensor networks; 
the operational resources required in the simulation of an 
n-qubit system scale like Poly(n, 2 X ) 16] (c.f. Theorem 
4 therein). Thus, MQC on L x L-planar code states, with 
X growing at least linearly in L, is not suited for efficient 
simulation by the methods developed in |3] - (HI- 

B. The 2D cluster state in the Ising model 

There is a two-dimensional universal cluster state as- 
sociated with an edge-centered graph (see Fig. whose 
overlap with an arbitrary product state is described by 
the Ising model partition function. However, its interac- 
tion graph is not planar. 

The cluster state \Ci) under consideration is shown in 
Fig. 01 a. Qubits are associated with the vertices of the 
displayed graph Cl- The stabilizer of Cl is generated by 
the Pauli operators 

K a = X a (g) Z b , VaeV{C L ). 

b\ (a,b)£E(C) 

The graph Cl is bi-colorable. We refer to the two colors 



as "i?" and "£?" , and use symbols "0" and V to dis- 
play the respective vertices. The cluster state \C T ) is a 
universal resource for quantum computation • 

Fig-El 3 shows a surface code state \Gl) based on non- 
planar graph Gl- Qubits are associated with edges in 
Gl- Gl is a modification of the previously discussed 
L x L-lattice graph. It has one extra vertex vo, extra 
edges (t>o, Vi) and extra plaqettes (vq, Vi, Vj), where Vi, Vj 
are vertices in the plane incident to a common edge. 
The stabilizer of |Gl) is spanned by the operators A s , 
Bp defined in Eq. QJ. For an L x L-lattice there are 
N{ = 2L(L + 1) + (L + l) 2 qubits, N[> = (L + l) 2 inde- 
pendent site- and N% = 2L(L + 1) independent plaquette 
stabilizers, and \Gl is thus uniquely specified. \Gl) may 
again be written in the form of Eq. © , and its overlap 
with a multi-local state may again, as in Eq. © be re- 
lated to the partition function of the Ising model. How- 
ever, the underlying interaction graph Gl is now non- 
planar such that the techniques for efficient simulation 
j26|,|27| are no longer applicable. 

The cluster state |Cl) and surface code state |Gl) are 
local unitary equivalent. First, they have support on 
the same number of qubits, and we can pairwise identify 
the respective qubit locations. The mapping is as fol- 
lows. The vertices of color R in C reappear as vertices of 
Gl, but they are no longer qubit locations. As for the 
latter, each i?-colored vertex v of Cl corresponds to an 
edge (v,vq) of Gl above the plane, and each i?-colored 
vertex v' of Cl, located between two i?-colored vertices 
Vi, V2, corresponds to the edge (vi,t>2) of Gl inside the 
plane; See Fig. 0] Denote by Hr the simultaneous local 
Hadamard transformation on all qubits of color R. Then, 
\G L ) = H R \<f>)c L - 

To see this, first consider the stabilizer K w of |Cl) 
associated with a -B-colored qubit location w, K w — 
(vi, v 2 are vertices of color R). Under the 
identification of vertices in Cl with edges in Gl, Wi 2 
{vo,vi t 2), w *-* (1)1,1)2). (Again, vq is the extra vertex 
above the code plane.) After application of Hr, the op- 
erator K w is mapped to X^ Vl , V2 )X^ Vo , Vl )X {vo , V2 ) = B p 
for p = (vo, Vi,V2); c.f. Eq. 10} . Second, consider the 
stabilizer K v for an i?-colored vertex v £ V(Cl), K v = 
X v ® w \ ( v ,w)£E(c) Z w . With the identification v <-> (v, vq), 
w <-> (v, w), and after application of Hr, the stabilizer K v 
is mapped onto Z {VtVo) ® w \ {v , w)eE (G L ) z {v,w) = A v . □ 

The overlap between a local state |$) = <^> eG £; \<f> e ) 
and |Gl) may again be written as the partition function 
of the Ising model, 

(Gl|$) - J2 CX P [H PjkVjCTk 

Yl ex P I H (3jk<Tj<Jk + J2/3joVj ■ 

(13) 



8 

v 




FIG. 5: Intermediate state between the planar code state and 
a two-dimensional cluster state. 

Therein, the (3jk are specified through the relation 

exp(fijk) = \J (0\4>jk)/(l\4>jk), and \<p e ) = H e \<j> e ). As 
can be seen from the lower line in Eq. (|13fl . the overlap 
between a local state and a universal 2D cluster state 
corresponds to the partition function of the planar Ising 
model in the presence of a magnetic field. 



C. The complexity of MQC simulation for 
non-planar interaction graphs 

First, we would like to morph the simulation of MQC 
with a planar code state into simulation of MQC with 
the universal cluster state, by 'switching on magnetic 
fields'. This may be performed gradually, switching on 
one magnetic field at a time. Fig. shows the graph 
corresponding to an intermediate state of this sequence. 
The sequence is passed in reverse order if one starts from 
a 2D cluster state of Fig.0]and measures, one by one, the 
cluster qubits of color R in the A-basis. More generally, 
X- and Z-measurements directly operate at the level of 
Ising interaction graphs G; by deletion and contraction 
of edges, respectively. 

A measure for deviation from a planar graph is the 
number rj of edges that needs to be deleted from a graph 
to make it planar. This parameter rj also governs the 
complexity of a straightforward extension of the simula- 
tion scheme presented in Section UTTI applicable to non- 
planar graphs. Be E np a set of edges in G such that 
G\E np is planar. Then, the product state &) eeE \4>e) 
may be expanded into the A-basis and the efficient sim- 
ulation of MQC on the remaining planar graph may be 
run for each component. The operational cost of this 
simulation method for an rt-qubit state |G) is Poly(n)2 v , 
which may be compared to Theorem 4 of [l6| for tree 
tensor networks. 



V. CONCLUSION 

We have shown that the overlap between the planar 
code state and any product state can be computed ef- 
ficiently classically, through its correspondence with the 
planar Ising model. Under rather general assumptions 
about the pattern of one-qubit measurements, MQC with 



9 



the planar code state as a quantum resource is not uni- 
versal. It can be efficiently simulated on a classical com- 
puter. MQC with a universal 2D cluster state can also be 
related to the Ising model. However, the corresponding- 
interaction graph is non-planar. 

Acknowledgments 

We are thankful to Alexei Kitaev for useful discussions 
and reading the manuscript. RR would also like to thank 



Panos Aliferis, Maarten van den Nest and Herbert Wag- 
ner for discussions. SB acknowledges support by the NSA 
and ARDA through ARO contract number W911NF-04- 
C-0098. RR is supported by the Government of Canada 
through NSERC and by the Province of Ontario through 
MEDT. 



[1] R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982); Opt. 
News 11, 11 (1985); Found. Phys. 16, 507 (1986). 

[2] G. Vidal, Phys. Rev. Lett. 93, 040502 (2004). 

[3] R. Cleve, A. Ekert, Ch. Macchiavello, M. Mosca, Pro- 
ceedings of the Royal Society of London, Series A 454, 
No. 1969, 339 (1998). 

[4] G. Giedke and J.I. Cirac, Phys. Rev. A 66, 032316 
(2002). 

[5] S. Bravyi, Quant. Inf. Comp. 5, Vol. 3, 216 (2005). 
[6] B.M. Terhal and DP. DiVincenzo, Phys. Rev. A 65, 
032325 (2002). 

[7] D. Gottesman, Proceedings of the XXII International 
Colloquium on Group Theoretical Methods in Physics, 32 
(1999). 

[8] R. Raussendorf and H.-J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[9] N. de Beaudrap, V. Danos, E. Kashefi, quant- 

ph/0603266. 
[10] W. Hall, quant-ph/0512130. 
[11] D. Gross and J. Eisert, quant-ph/0609149. 
[12] A.M. Childs, D.W. Leung, M.A. Nielsen, Phys. Rev. A 

71, 032318 (2005). 
[13] M. van den Nest, A. Miyake, W. Diir and H.-J. Briegel, 

Phys. Rev. Lett 97, 150504 (2006). 
[14] I. Markov, Y.-Y. Shi, quant-ph/0511069. 
[15] Y.-Y. Shi, L.-M. Duan and G. Vidal, quant-ph/0511070. 
[16] M. Van den Nest, W. Diir, G. Vidal, and H. J. Briegel, 

quant-ph/0608060. 
[17] M.A. Nielsen, quant-ph/0504097. 
[18] A. Kitaev, Ann Phys. (NY.) 303, 2 (2003). 
[19] S. Bravyi and A. Kitaev, quant-ph/9811052. 
[20] R. Raussendorf, S. Bravyi and J. Harrington, Phys. Rev. 

A 71, 062313 (2005). 
[21] E. Dennis, A. Kitaev, A. Landahl, J. Preskill, J. Math. 



Phys. 43, 4452 (2002). 
[22] A. Hamma, R. Ionicioiu, and P. Zanardi, Phys. Lett. A 
337, 22 (2005). 

[23] S. Bravyi, M. B. Hastings, and F. Verstraete, Phys. Rev. 

Lett. 97, 050401 (2006). 
[24] A. Kitaev, private communication (2005). 
[25] M. van den Nest, W. Diir and H.-J. Briegel, quant- 

ph/0610157. 

[26] F. Barahona, J. Phys. A: Math. Gen. 15, 3241 (1982). 
[27] P. Kasteleyn, Physica 27, 1209 (1961). 
[28] A. Galluccio and M. Loebl, Electron. J. Combinat. 6, R6 
(1999). 

[29] J.K. Pachos, M.B. Plenio, Phys. Rev. Lett. 93, 056402 
(2004). 

[30] D. Schlingemann, R.F. Werner, Phys. Rev. A 65, 012308 
(2002). 

[31] M. Grassl, A. Klappenecker and M. Rotteler, IEEE in- 
ternational symposium on information theory, Lausanne 
(2001). 

[32] M. Hein et al, Proceedings of the International School 
of Physics "Enrico Fermi" on "Quantum Computers, 
Algorithms and Chaos", Varenna, Italy (2005); quant- 
ph/0602096. 

[33] S.-I. Oum, PhD-thesis, Princeton University (2005). 

[34] Note that at each step of MQC only some part of qubits is 
measured. To compute probabilities of outcomes of such 
partial measurements we have to consider arbitrary sub- 
graphs of the square lattice. 

[35] One can easily check that projecting a qubit onto a state 
0) or |1) effectively removes this qubit from the lattice 
leading to a state \Ge) on a modified planar graph. One 
can use this observation to avoid the cases a e — 0. 



